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We propose a protocol for measurement of the phonon number distribution of a harmonic oscillator 
based on selective mapping to a discrete spin-1/2 degree of freedom. We consider a system of a har¬ 
monically trapped ion, where a transition between two long lived states can be driven with resolved 
motional sidebands. The required unitary transforms are generated by amplitude-modulated poly¬ 
chromatic radiation helds, where the time-domain ramps are obtained from numerical optimization 
by application of the Chopped RAndom Basis (CRAB) algorithm. We provide a detailed analysis 
of the scaling behavior of the attainable fidelities and required times for the mapping transform 
with respect to the size of the Hilbert space. As one application we show how the mapping can be 
employed as a building block for experiments which require measurement of the work distribution 
of a quantum process. 

PACS numbers: 02.60.Pn,03.65.Aa,03.65.Wj,32.80.Qk 


I. INTRODUCTION 


Trapped ions represent a system where continuous and 
discrete degrees of freedom can be jointly manipulated 
and measured in the quantum regime m- Continu¬ 
ous degrees of freedom are given by oscillatory motion 
in the trap potential, while the discrete ones are given 
by internal (spin) states. For coherent spin manipula¬ 
tions, optical or microwave radiation can be employed. 
This enables a wide variety of applications in the fields 
of quantum computing ma, quantum simulation PIT] 
and quantum metrology [5]. In this work we propose a 
scheme that enables single-shot interrogation of the har¬ 
monic oscillator degree of freedom and thus paves the 
way for novel applications of trapped ions. This is mo¬ 
tivated by the persisting difficulty to implement single¬ 
shot and/or quantum-non-demolition (QND) measure¬ 
ments of the phonon number for trapped ions, which is 
mainly due to the predominant harmonicity of the trap 
potential. Ref. P proposes a QND filtering scheme, 
while Ref. [10] demonstrates a single-shot readout. Both 
schemes require multiple iterations of unitary manipu¬ 
lation and spin readout, which is time consuming and 
imposes significant experimental complexity. While we 
consider the paradigmatic system of trapped laser-cooled 
ions, the scheme is applicable to other systems such as 
atoms in cavities midi] or superconducting qubits m- 
It is therefore of general interest to study the quantum 
controllability of this class of systems in detail and under 
the assumption of realistic parameters. 

In this work, we exploit the controllability of the 
system to construct a selective, unitary phonon-to-spin 
mapping scheme, where the spin degree of freedom un¬ 
dergoes a flip operation conditioned on the phonon num¬ 
ber of the motional mode. Our scheme is fully unitary, 
however it does not directly serve as a single-shot phonon 
number measurement, it rather yields a dichotomic result 


- namely whether the system is found in a specified num¬ 
ber state or not. This is a consequence of the fact that the 
spin-1/2 degree-of-freedom is probed. We show how this 
can be enhanced to yield a QND measurement scheme. 

The controllability of the problem has been analyzed 
by [II], and the control problem has been tackled for 
state-to-state transfer m and gate optimization [TSIITB]. 
This work is an example for incomplete control: We seek 



FIG. 1. Schematic depiction of the spin mapping Um for the 
case of m = 2 (initial state left, final state right): Starting 
from a product state with spin down and occu¬ 

pation in different number states of the harmonic oscillator 
(the circle diagrams show |c„p), only the occupation in the 
number state |n = m = 2) (|c 2 p, yellow) shall be excited to 
lb) while the occupation initially in number states \n ^ 2) 
knP, gray) shall end up in |/). As shown on the right 
after the mapping |c 2 p (yellow) is distributed over the spin 
up states while (gray) is distributed over the spin 

down states. Note that only the spin population of the fi¬ 
nal state matters while the state of the motional mode is not 
specified. 

to map a specific vibrational level m S 0... fV — 1 on 
a spin state, see Fig. In other words, we construct 
a unitary transformation which flips the spin only for 
a predefined vibrational number state m and leaves the 
spin unaffected for the remaining levels n ^ m. However, 
the state of the motional degree of freedom for the final 
state, is not specified. This allows us to fix only N pa¬ 
rameters (specifying the final time spins conditional the 
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FIG. 2. Performance of Um{T) for different m as it acts on the 
states I 4,,n). a) m = 0, b) m = 1, c) m = 2, d) m = 3. The 
columns show the resulting population in the upper spin state 
p-f = |(t, k\Um{T)\n,X)\^ depending on the initial phonon 

number state |n): the yellow colums show the population for 
m, where the spin should be flipped, the gray columns show 
the spin up population for phonon numbers n, where we want 
to keep the state in spin down. 


initial harmonic oscillator state) of the full 4fV^ — 1 di¬ 
mensional unitary, thus reducing the complexity of the 
control problem and allowing for higher fidelities in a 
shorter time. In this sense the control is incomplete and 
we can treat systems with larger N. 


Already for small sizes of the respective Hilbert space 
no intuitive solution for this control problem is available. 
We therefore construct the required control fields numer¬ 
ically, based on the Chopped Random Basis (CRAB) 
[niiiH] optimization algorithm. 


In this work, we first outline the mapping scheme, and 
integrate this mapping into a quantum non demolition 
filter in section |lll In section lTTTl we will describe in detail 
the system of the trapped ion that we want to use to real¬ 
ize the filter, and define the fidelity measure of the map¬ 
ping. 


Section IV will introduce the numerical method 


and parameters for simulation of the system and section 
[V| will specify the optimization method and present re¬ 
sults from optimization. 


A. The Mapping Scheme 


The purpose is to map a specific vibrational level \m) 
on a spin degree of freedom, see Fig. We start with 
a product state p(0) = 1),) ® p^®(0), where the spin is 
initialized to j),) and the harmonic oscillator mode is in 
an arbitrary pure or mixed state characterized by (0). 
By applying a pulse sequence of up to three radiation 
fields we create a unitary transformation Um, that flips 
the spin for initial occupation in the vibrational level \m) 
and keeps spin down for |n) ^ \m). That is, Um maps 


ID® 




n = m 
n ^ m 


( 1 ) 


We will distinguish the ideal map Um that we seek to 
implement from the actual finite-time evolution Um{T). 
Fig.@ shows the sucess of such maps for m = 0,1, 2, 3,4 
as obtained by CRAB. The details of the results will be 
explained in section [V] 

Note that for each initial state |n), the vibrational 
mode is supposed to be in an unspecified superposition 
state after the filter unitary transform. Thus, only direct 
readout or coherent spin manipulations are to be applied 
after the filter operation. Using trapped cold ions, the 
spin can be read out with fidelities of 99.99% [TS], and 
preparation of identical input states is current state of the 
art. We can repeat a sequence of preparation, mapping 
transform Um(T) and spin readout to obtain a statistical 
estimate for pm = P^m^ the initial occupation in |m). 
Performing this for all number states in a truncated sub¬ 
space yields the phonon distribution. With additional 
displacement operations after preparation, this can be 
extented to a full tomography scheme m- 


II. APPLICATION FOR FILTERING AND THE 
MEASUREMENT OF WORK DISTRIBUTIONS 

In this section, we outline how the unitary mapping 
operation can be used in conjunction with spectroscopic 
decoupling [2T] to realize a filtering protocol similar to 
the QND filter proposed in 0 US- In that scheme, 
pulses of defined area on a motional sideband realize 
a low-precision mapping operation. After one mapping 
operation, the spin is projectively measured by fluores¬ 
cence detection. Absence of fluorescence indicates that 
the system is likely to be in |m), and the discrimination 
efficiency increases upon iterative repetition. This corre¬ 
sponds to preparation in level |m) by QND measurement. 
For a quantum process acting on the vibrational mode, 
this scheme allows for measurement of the work distri¬ 
bution by QND preparation of level |m), the subsequent 
quantum process C and a second filtering operation for 
level \m'). This yields the probability that the quantum 
process transfers level |m) to |m'). 

While this opens up the prospect of studying the non¬ 
equilibrium thermodynamics of a well-controlled quan¬ 
tum system, the approach is difficult to implement ex¬ 
perimentally. The main reason lies in the fact that a 
Fock state needs to be preserved over a comparatively 
long time span, where the repeated unitary driving and 
dissipative measurement operations take place. The fil¬ 
tering scheme presented in Ref. uni works by iterative 
adiabatic phonon removal, where the event that the sys¬ 
tem has initially been in Fock state |n) is heralded by 
bright detection after n — 1 previous cycles of phonon re¬ 
moval and dark detection. The scheme however requires 
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FIG. 3. Employing unitary phonon to spin mapping for a QND filter: The panels show a sequence of mapping and spectroscopic 
decoupling operations which enables measurement of the work distribution of a quantum process, a) shows a state prepared 
in 14,) with an arbitrary state of the vibrational mode, b) shows how population from |m = 2) is transferred to |t) by means 
of C/ 2 . The vibrational populations are reshuffled in an unspecihed way. In c), population from |4,) is shelved to a metastable 
auxiliary state |ai) for spectroscopic decoupling, d) shows how the population which was not shelved is transferred to |4,, 2) 
by applying the time reversed mapping C/J. In e), a general quantum process £ (unitary and/or dissipative) takes place which 
performs work on or transfers heat to the vibrational mode. In f), the population of \m' = 0) is mapped to |t) by means of 
Uo, and the remaining population in |4,) is shelved to another metastable state |a 2 ) in g). Finally, h) shows how readout takes 
place by cycling to population of the ground state through an excited state |e) and detecting the resonance fluorescence. 


the deterministic preparation of state \n) after the final 
bright event, e.g. by convential ladder climbing schemes 
[23l[2l]. By contrast, our scheme consists of a fully uni¬ 
tary mapping operation, which can be enhanced to a 
measurement protocol for work distributions by popu¬ 
lation transfer from a specific spin state to additional 
met a-stable levels. 

The scheme is explained in detail by Fig. We start 
out with an initial state |4.) ® /O^°(0). We might be in¬ 
terested in the population of some vibrational level \m) 
(panel a), m = 2). We apply the mapping Um to transfer 
this population to |t) (b), and the population remaining 
in 14.), i.e. for all |n y/ m), is transferred to the metastable 
state |ai) (c). We then apply the time reversed filter 
and obtain the state |4,, m), ignoring the shelved popula¬ 
tion (d). We then carry out the generic quantum process 
C which is to be analyzed, which acts only on the mo¬ 
tional degree of freedom (e). The redistribution of popu¬ 
lation among the Fock states is then analyzed by applying 
the mapping operation Um' (f) and shelving the popula¬ 
tion left over in |4,) to a second metastable state [ 02 ) 
(g). We can now measure the groundstate population by 
detection of state-dependent resonance fluorescence by 
driving a cycling transition to a short-lived excited state 
|e) (h). For a quantum process p' = C{p), the measure¬ 


ment sequence yields fluorescence with probability 

Pf = (m|/9^°(0)|TO) (m'|£(|m)(TO|)|m'). (2) 

The initial populations can be measured separately by 
simply performing the fluorescence detection after the 
first shelving operation. Thus, the measurement scheme 
yields the transfer matrix of the quantum process C and 
therefore its work distribution. 

The scheme can be implemented with all ion species com¬ 
monly used in experiments. For or ®®Sr+, the 

states It), It) can be identified with the Zeeman sublevels 
of the Si /2 ground state, the auxiliary states |ai),|a 2 ) 
with different sublevels of the 1 / 5/2 state, which has a 
lifetime of about 1 ms, while the excited state would be 
the Pi /2 state. 


III. THE SYSTEM 

We consider a trapped ion of mass M and limit our¬ 
selves to a ID harmonic oscillator of frequency The 
internal states are coupled by radiation, which is char¬ 
acterized by the Rabi frequency fl. A sufficiently short 
wavelength A serves to couple internal and vibrational 
levels. This condition is described by the Lamb-Dicke pa- 
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tional (harmonic oscillator) degrees of freedom. Up to 
three radiation fields are applied to drive transitions be¬ 
tween these levels: One resonant with the carrier transi¬ 
tion which connects the two spin states for equal phonon 
number, and fields driving the blue (red) sidebands. The 
combined phonon and spin level scheme as well as the 
laser driven transitions and the dependence of the rela¬ 
tive transition strength on the vibrational excitation n 
are depicted in Fig. for the parameters given in the 
following section. 


FIG. 4. a) Energy level scheme of a single trapped ion. The 
relevant levels are product states of the spin state and the 
harmonic oscillator modes given by the phonon number states 
|n). Lasers drive transitions on the carrier as well as on the 
blue and red sidebands, b) Relative Rabi frequencies (given 
by the matrix elements Mn,n+Sn = (^)) 

T] — 0.25. On resonance this is the relative coupling strength 
of the carrier transition (black circles, Sn = 0), the blue (stars) 
one for the blue sideband {5n = 1), and the red (crosses) one 
for the red sideband {Sn = —1). 


rameter V = jy 2 Muj ’ specifying the ratio of the ground 
state wavepacket size to the wavelength of the driving ra¬ 
diation [25] . 

Any ion species which features a closed transition for 
laser cooling and a qubit transition is suitable. Thus, op¬ 
tical m, stimulated Raman [26] and microwave qubits 
m can be employed. The Hilbert space of the system is 
a product space of the two relevant spin and the vibra¬ 


A. Hamiltonian 

Three radiation fields contribute to the Hamiltonian. 
We label the fields by a = c, 6 , r and characterize them by 
the bare Rabi frequency and the detuning from the 
carrier transition Sc°‘'^ = uja—u!c {tOc the carrier transition 
and oJa the frequency of the field). To have the fields on 
resonance with the carrier (Hq'^^) and the blue and red 
sidebands (Hg**^ and Hg’^^ respectively) we have to choose 
= 0 , Si^'^ = and = —oj^, where the trap 
frequency ujz is the energy of one phononic excitation. 

Each radiation field does not only resonantly drive 
transitions, it also causes energy level shifts via off- 
resonant driving. We neglect terms rotating at 2uJz or 
faster, and terms which scale with a |5n| > 1 matrix 
element. The Hamiltonian, including resonant and off- 
resonant effects, then reads: 




n=0 


(f\ 

-> + fa)(n| 

+ 5n) {n\ 

ol'') (t) 


-I- h.c. 


( 3 ) 


The Rabi frequencies depend on the phonon number n 
and the Lamb-Dicke factor 77 as 

= (4) 

where the matrix elements are modified Laguerre poly¬ 
nomials and are given in the appendix. This dependence 
of the Rabi frequency on the phonon number is the cor¬ 
nerstone for any conceivable filter mechanism. 

Example matrix elements are shown in Fig.|^b). Note 
that the differences between the matrix elements become 
smaller for increasing n, thus larger discrimination times 
are required as the relevant part of the Hilbert space (i.e. 
if we want to have our map acting in the desired way for 
a higher range of initial states) increases. 


B. Fidelity Function 

Since the ion will initially be in a state with only a few 
motional excitations, we evaluate the fidelity of Um{T) 
(the imperfect implementation of Um) only by its ac¬ 
tion on an Wdimensional subspace spanned by the basis 
states {|7/>n)}n=o...Af-i, with |■^/'„(0)) = | i,n). The map 
Um{T) evolves |V'n(0)) into \il’n{T)) = 17^(T)jV'„(0)). We 
choose the fidelity F of Um{T) to be the product of the 
probability that the desired spin flip occurs for level to, 
and the averaged probability that no undesired spinflip 
occurs for n ^ m: 


Fim,N)=F^F^, 


( 5 ) 
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where 


N-l 

F^=Y,\{t,k\^^iT))\\ ( 6 ) 

k^O 

N-l 

^^ = ]vr3i E E (7) 

n^m k—0 

So if = 1 it means that for any input state the spin up 
population after the mapping corresponds exactly to the 
population of \m) before the mapping. In other words 
the spin is fully flipped to |t) for initial state \m) and 
stays in jj,) otherwise. 

The maps Um will by design of F{m,N) produce the 
desired spin operation only for those initial states that 
do not populate vibrational levels \n > N). Instead Um 
will act in an unspecified way on |n > N). Thus N 
initial basis states are included when calculating the fi¬ 
delity. Throughout the article we choose N = 10, unless 
explicitely stated. We will then write to refer to the 
number of initial states included in the fidelity of the 
map. 


IV. SIMULATION AND PARAMETERS 




T(ps) 

FIG. 5. Pulse shapes after optimization. Pulse height in 
units of SIo = Stt • 50 kHz. a) The optimal pulse for m = 1, 
T — 400 /iS, where three radiation fields (black=carrier, 
red=red sideband, blue=blue sideband) are applied and all 
of them are modulated continuously in time. The transfor¬ 
mation works on the harmonic oscillator levels n = 0,... 9. 
b) Two fields (black=carrier, red=red sideband) are applied 
and the control is limited to phase switching by tt while the 
power is kept constant = S7o or = —Ho)- 

The transformation operates on the harmonic oscillator levels 
n = 0,... 3 for m = 0, T = 300 fis. 


The simulation of the unitary dynamics is carried out 
by truncating the Hilbert space, allowing for vibrational 
excitations of up to \n = 14), i. e. the total dimension is 
30. The Hamiltonian is mapped on this subspace and 
time dynamics can be simulated by diagonalizing the 
Hamiltonian at each time step. The validity of the trun¬ 
cation is verified by checking convergence of the simula¬ 
tion on a bigger subspace with dimension 40. 

Realistic parameters for the ion light interaction are 
chosen with a maximal bare Rabi frequency Hq = 27r • 
50 kHz, operation time T = 100— 1000 /rs with 1000 time 
steps that determine the grid for changes of the field am¬ 
plitudes and thus allow us to represent pulse modulation 
frequencies up to about (10 — 20) • 27r/r. The Lamb- 
Dicke parameter was set to 77 = 0.25, the trap frequency 
Wz = 27r • 1.4 MHz, so well in the regime ujz 3> Hg where 
resonant terms are dominant over offresonant terms. The 
parameters correspond to recent experiments [231126] . 


V. OPTIMIZATION 

Now we are ready to engineer Um by simulating the 
system’s time evolution and optimize the pulses via the 
Chopped Random Basis (CRAB) algorithm [TTHTH]. The 
Rabi frequency ramps lead to the dimensionless pulses 

/(“)(t) = H(“)(<)/Ho, (8) 

which are optimized so that the pulses vanish at zero 
and final time and |/*'“Ht)| < 1 , i.e. the maximum pulse 
height does not exceed the experimental constraint Og. 


A. Control with three radiation fields 

If we control the laser power on the carrier transition as 
well as on the red and blue sideband we have three (real) 
controls that are subject to CRAB. If we allow also for 
the phases of the lasers to be controlled by CRAB we 
end up with 6 (real) controls, but it turns out that this 
does not lead to substantial improvement of the fidelities. 
Therefore in the following we restrict ourselves to real 
valued pulses. Note that we assume the lasers to be in 
phase at t = 0 , although an initial phase different will 
not cause a big effect [28] . 

Optimization was done for selecting m = 0,1,2,3. 
Fig. shows how the achieved mapping unitaries UmiT) 
transform initial spin down population into spin up pop¬ 
ulation depending on the initial state |'0„(O)) along with 
the fidelities of the operations. 

In an actual experimental realization, the main error 
source would result from an imperfect calibration of the 
Rabi frequencies. This error can be modelled as 

fio/(“)(t) = Hg(l+0/^“^(t) (9) 

where Hg/(“)(t) are the actual Rabi frequencies while 
{t) is the calculated optimal pulse. The robustness 
of the results is thus characterized by the fidelity as a 
function of For —0.01 < ^ < 0.01 the fidelity does not 
differ by more than 1.3% from its maximum value. 
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B. Control with two radiation fields 

If we use only two driving fields the control task is 
more difficult to achieve. However, this is a very inter¬ 
esting problem not only because it requires less complex¬ 
ity when implemented in an experiment, but also since 
it corresponds to the systems analyzed e.g. by nnni- 
[Tfi] , Figure [^compares the resulting fidelity F{m, 10) for 


0.9 

0.85 


FIG. 6. Optimization results for different control scenarios: 
fidelity F(m, 10) for three driving fields (black squares), driv¬ 
ing only the carrier and the blue sideband transition (blue 
circles) and using only the carrier and the red sideband tran¬ 
sition (red crosses) as a function of m. 

three driving fields (black squares), with scenarios using 
the carrier and the blue sideband transition (blue circles) 
and using the carrier and the red sideband transition (red 
crosses). Both options work equally well. In the follow¬ 
ing we consider one field resonant to the carrier and one 
field resonant to the red sideband transition since it is 
closer to the systems in Ref. [nmiHis]. 

So far we had Um operating in the desired way on 
the subspace of = 10, thus for up to 9 vibrational 
quanta. This was achieved by maximizing the control 
functional F{m, N = 10) that considers the time evolu¬ 
tion of iV = 10 basis functions. However, in a scenario 
with only few excitations it can be enough to distinguish 
e.g. \n = 0) from \n = 1) and |n = 2) while it is known 
that there is no contribution from \n > 3). This requires 
a map that acts in a desired way only on |V'n)n=o,i ,2 and 
can be achieved by maximizing F(m,N = 3), thus con¬ 
sidering the time evolution of only iV = 3 basis functions. 
As presented in m, the control task dramatically simpli¬ 
fies if fewer levels are included in the control functional. 
Fig.j^shows the scaling of the fidelity F(m, N) with N for 
m = 0 and T = 300/rs (black squares). As one can see, 
especially for small N the fidelity is much higher than for 
intermediate and larger N. One obvious reason, as men¬ 
tioned already, is that we have to fix fewer parameters 
for smaller N (e.g. the time evolution of 3 basis functions 
for as compared to 10 basis functions in the case of 
U}^) while the complexity of the optimization problems 
scales with the number of parameters [29] . Furthermore, 
for larger N the transitions with higher n are included, 
where the matrix elements Mn^n+Sn of the carrier and 
the sidebands become comparable to one another and 
the scaling with n becomes smaller and thus distinguish¬ 
ing between two vibrational levels becomes harder. The 
influence of this scaling of the transition strength was 
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also studied in appendix [C] by considering the system as 
a set of Poincare pointers. 

A control strategy that is even simpler to realize ex¬ 
perimentally but also more limited would be to keep the 
pulses on constant power and modulate the pulses just 
by phase flips (H(“)(t) = Oq or 0(“)(t) = —Hq)- We call 
this a discrete pulse. The corresponding pulse is shown 
in Fig. [^b), while the resulting fidelities F{m = 0,N) 
are shown by gray circles in Fig. [7] We present results 
for A^ = 2,... 5 since the restrictions on the pulse make 
the control task much more challenging and if we include 
more levels the fidelity drops. Still we get good fidelities 
for small N. 

10-1 
l-F 

10-2 

10 '^ 

10-'' 

10 -^ 


FIG. 7. Operation error 1 — F as a function of N, the num¬ 
ber of basis states considered. The results are for realizing 
Uo^ in an operation time of T = 300 /rs using carrier and red 
sideband. The black squares show the results allowing con¬ 
tinuous modulation of the pulse power over time. The gray 
circles show the results obtained when only phase flips at con¬ 
stant power are allowed (compare Fig.for the pulse shape). 


Uq, continuous pulse ■ 

Uq, discrete pulse • 
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VI. CONCLUSIONS 

We have obtained control fields for phonon-selective 
flipping of the internal state of a trapped ion by appli¬ 
cation of the CRAB algorithm. The optimization yields 
solutions with fidelities of better than 95 % where total 
times < 500/is are required for the mapping transform 
laser pulse. This bears interesting applications for the 
quantum state reconstruction of the motional mode, and 
for the measuring of the work distribution pertaining to 
quantum processes. This in turn opens up perspectives in 
the emerging field of quantum thermodynamics [13130]. 
An important feature is that we do not specify the final 
state of the harmonic oscillator degree of freedom. This 
is an algebraic reduction of the degrees of freedom to be 
fixed by optimal control and has been shown to lead to 
significantly better control performance compared to full 
gate optimization in a slightly different context |31| . An¬ 
other interesting aspect of the mapping transform is that 
for only one laser it reduces to the Poincare recurrence 
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time problem. This analogy is studied in appendix 
While we focus on trapped ions the results are also rele¬ 
vant for systems of atoms in cavities and superconducting 
qubits. 


The fact that the Rabi frequencies depend on the vi¬ 
brational quantum number n enables our desired spin 
flipping mechanism. 
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Appendix A: Hamiltonian 


First we consider the contribution of a single radiation 
field to the Hamiltonian. We label the fields hy a = c,b,r 
and characterize them by the bare Rabi frequency 
and the detuning from the carrier transition = Wq, — 
LOc (wc the carrier transition and the frequency of the 
field). To have the fields on resonance with the carrier 
and the blue and red sidebands and 

respectively) we have to choose = 0 , = ujz, and 

= —uiz, where the trap frequency uJz is the energy of 
one phononic excitation. 

The contribution of one of these fields to the total 
Hamiltonian in the interaction picture reads 




n^O ^ 

+ l)(n| + h.c. 

(Al) 


thereby neglecting contributions that connect phonon 
numbers that differ by more than one. The Rabi frequen¬ 
cies depend on the phonon number n and the Lamb-Dicke 
factor T] as 

= (A2) 

where the matrix elements are 

MnAv) = (77^) (A3) 

Mn,n+i{v) = (A4) 

V n + 1 

Mn,u-i{ri) = . (A5) 

The functions are the standard Laguerre polynomials. 

We assume that the amplitudes of the drive fields (t) 

can be varied arbitrarily over time. 

The total Hamiltonian then reads 

H{t) = (t) + (t) -f (t). (A 6 ) 


Appendix B: OfFresonant terms — Stark shift 


As we can see in equation (Al) each laser contributes 
not only the resonant term but also non-resonant terms. 
We neglect terms rotating at 2w^, and terms which scale 
with a |<5n| > 1 matrix element. The remaining four 
terms are: 


^ n,n±l -| ^ |n ± 1) (n| F.C. , (Bl) 

n—0 

Clip) 

I 0 |n)(n| -I- h.c. and(B2) 

n—0 

E t)(; I ® \n){n\ + h.c. ,(B3) 

n—0 


Together with the resonant terms this corresponds to 
equation All other offresonant terms are neglected 
and even these four terms are small since uJz ^ 


Appendix C: The Poincare Recurrence Time 

The desired spin mapping operation Um could be ac¬ 
complished also without sophisticated control pulses but 
with a constant pulse on the e.g. blue motional side¬ 
band. This is related to Poincare’s [35] recurrence time 
theorem. The theorem states that a conservative, finite, 
closed system returns arbitrarily close to its initial state 
after the Poincare recurrence time. Hemmer et al. |33j 
showed that for a chain of N harmonic oscillators this 
boils down to the problem of N pointers, all rotating 
with a certain frequency a;„, where these frequencies are 
incommensurate. Thus, they can calculate the Poincare 
recurrence time from the probability of having all point¬ 
ers directing in the same direction up to a directional 
error A^jj. 

If we try to solve our spin mapping problem for the 
simple case of driving the blue sideband (where we can 
choose the amplitude to be constant without loss of gen¬ 
erality) on the subset of n = 0 ,.., A — 1 , this corresponds 
to a set of N rotating pointers, that in the end all have 
to point in a certain given direction (given by the desired 
spin flip described by Um) up to an error A(^„. The fre¬ 
quencies of the pointers are ujn = In the Lamb- 

Dicke regime rj'^n <C 1, ^n'^n+i ^ \/n + l, from which 
we see that the frequencies are actually incommensurate. 
Following [33] , the theoretical waiting time Tp after which 
Um is generated up to some specified error is 

AT-l N 

s*™+. n ^ (Cl) 

^ ^ ' n=0 fe = 0 

k^n 
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FIG. 8. Waiting times in units of Theoretical 

Poincare recurrence times Tp (red crosses) and required wait¬ 
ing times found by numerical simulations (black crosses and 
blue squares) are plotted for different numbers N of in¬ 
cluded phonon excitations. The results were obtained for 


= 27r • 50 kHz and e = 0.02. Blue squares are wait¬ 
ing times Tp for crossing the error threshold (equation |C5[ ), 
black crosses are waiting times for constant Acfi (equation 


C6l. The gray circles show Topt, the operation time needed 


for the optimized maps Uo{N) with carrier and red sideband 
presented in Fig. 


The errors due to A(j)k 7^ 0 correspond to a fidelity of 


Tt = cos^ A(j)^ 

n^m 

n^m 


(C2) 

(C3) 

(C4) 


where the approximation holds for ^ 1. For a given 
maximal error e = 1 — F thus we have to find A(^„ such 


that 


+ E 


n^m 


N-l 


< e, 


while under the further constraint of A(/)„ not depending 
on n the condition on 1 — F" translates into 


A(j)n < a/ e/2 n = 0 ,... A — 1 . 


(C 6 ) 


We calculate Tp with this latter condition (i.e. setting 
A(l)n = 


1 


/eAAyi ^-1 

o 


(27r) 


(C7) 


n—0 


where we set the operation error to £ = 0.02. 

This theoretical Poincare recurrence time can be com¬ 
pared to waiting times obtained by numerical time evo¬ 
lution: we evolve the system under constant (t) = 
2 tt ■ 50 kHz and stop when our conditions on the phase 
deviations or fidelity holds. We call Tp the time p oint , 
where in the numerical time evolution equation (C5) 


holds for the first time and likewise Ta^ the time point 
where equation (C6) holds for the first time. Figure!^ 


shows Tp (red crosses) together with Tp (blue square^ 
and Ta 0 (black crosses) for A = 2,... 5. For the nu¬ 
merical values different points for same A correspond to 
m = 0,... A — 1. The Poincare waiting time is an ap¬ 
proximation for the stricter condition of equation ( |C6[ ) 
and thus closer to Ta^. As seen in the figure indeed 
the condition on A(/)„ is stricter, meaning longer waiting 
times {Ta 4 , > Tp). 

These waiting times are now compared to the control 
scenario of two radiation fields (carrier plus red side¬ 
band) and the pulse shapes are optimized for each A 
as in Fig. The gray circles show Topt, the operation 
time required to achieve Uo{N) up to an error £ = 0.02. 
While the presented values are too few to make general 
statements about the scaling of T with A it is clear that 
by optimization one can gain a speed-up of several orders 
of magnitude. 
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